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We consider here chiral symmetry breaking through nontrivial vacuum structure 
with quark antiquark condensates. We then relate the condensate function to the 
wave function of pion as a Goldstone mode. This simultaneously yields the pion also 
as a quark antiquark bound state as a localised zero mode in vacuum. We illustrate 
the above with Nambu Jona-Lasinio model to calculate different pionic properties in 
terms of the vacuum structure for breaking of exact or approximate chiral symmetry, 
as well as the condensate fluctuations giving rise to a mesons. 
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I. INTRODUCTION 



Nambu and Jona-Lasinio (NJL) had linked chiral symmetry breaking |l| with properties 
of hadrons quite sometime back with pion as the Goldstone mode However, pion is both 
a quark antiquark bound state and a Goldstone mode. Hence through Goldstone theorem 
pion state along with its wave function as a quark antiquark pair should also get related 
to the vacuum structure. It is surprising that this aspect is absent in the very extensive 
literature on the topic 

We consider phase transition as a vacuum realignment with an explicit structure and use 
the techniques developed earlier [||. In this manner for potential models we had obtained 
the gap equation to be the same as derived through Schwinger Dyson equation. Through 
Goldstone theorem [0] here we had the extra feature that the pion state as a space localised 
quark antiquark zero mode of destabilised vacuum also gets determined @J^. We further 
discussed the effects of approximate symmetry breaking where the gap equation changes 
giving rise to a change in the pion wave function. We discuss the same here for Nambu 
Jona-Lasinio model. The reason for doing the same in NJL model is its mathematical 
simplicity and its present relevance in the context of Salam Weinberg symmetry breaking 
and top quark mass P,p!0|. 



We organise the paper as follows. In section II we consider the vacuum structure with 
quark antiquark pairs using an ansatz for the same by minimising the energy density. This 
gives rise to the conventional gap equation and involves a new description of phase transition 
with an explicit construct for the destabilised vacuum |^ . In section HI we identify the pion 
as Goldstone mode and relate its wave function with functions associated with the vacuum 
structure. In section IV we consider the vacuum structure for NJL model when chiral 
symmetry is approximate. We also derive here some familiar results of current algebra in 
the present framework. In section V we calculate the pion charge radius using the wave 
function determined from the vacuum structure. In section VI we consider fluctuation of 
the condensate mode to give a qualitative identification of a meson. Section VII consists of 
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discussions. 

The method here consists of using equal time algebra 
the ground state through a variational principle P,|6|,p^ . 



I3| along with construction of 



II. CHIRAL SYMMETRY BREAKING AND VACUUM REALIGNMENT 



We shall now proceed in the same manner as earlier |^ for the NJL model. Let us start 
with the effective Hamiltonian 



nix) = i^ixfi-ta ■ ^MxY + /■rf#;^(f)^^(f)V:^;;',(f - y)V^^(y)V5(y), (1) 



which has chiral invariance. In the above i,j stand for color indices and a, (3 stand for the 

ij,kl I 



spinor indices and V^^J^L{x — y) is the potential. For effective QCD based vector potential 



we may take 



(2a) 



where A" are the Gellman matrices. We may also have NJL model when we take the contact 
potential as 



-y) = G [{l'U{l\s6''6'' - {l'l'U{l'l\s{r'^UT'^)u\ X 6{x - y). (2b) 

Here G is the dimensional interaction coupling constant and 's are the isospin matrices. 
The field operators ip{x) may be expanded as 
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iIj{x) 



(27r)3/2 

where U and V are given by 



e'-^'^dk 



(3) 



Ur{k) 



ya-kj 



Uir] Vs{-k) 



-o 



Vis 



(4) 



for free chiral fields. The perturbative vacuum is defined by this basis when we have 
c/ I vac >= 0= c\ I vac >. We next consider a trial vacuum state given as 
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vac' >= U I vac >= exp{B^ — B) \ vac > 



(5a) 



with 



= J cir{kf uir'^ {a.k)visCisi—k) f{k)dk. 



(5b) 



Here f{k) is a trial function associated as above with quark anti-quark condensate. We may 
recall a similar construction in Bogoliubov Valatin approach We shall minimise the 

energy density of | vac' > to analyse the possibility of phase transition from \vac > to 
\vac' >. For this purpose we first note that with the above transformation the operators 
which annihilate | vac' > are given as 



bj{k) = Uci{k)U-\ 



(6) 



which with an explicit calculation yields the Bogoliubov transformation 



^ birik) ^ ^ 



Jls\ 



-k) 



cosf 



y jjisinf{a 



cosf 



J 



Clr{k) 

\cis{-k) j 



(7) 



Here a^-s = ujr^^a ■ k)vis. Using the above transformation (6) or (7) the expectation value 
of the Hamiltonian with respect to | vac' > is given as 



S =< vac I H{x) I vac' >= T + V, 



(8) 



where T and V are the expectation values corresponding to the kinetic and the potential 
terms in Eq.(l). With a straightforward evaluation we then obtain that 

2N 



T =< vac' I tjj''{xy{-ia.\/)4)\x) \ vac' >- 



dk \ k \ cos2f{k), 



(9) 



where N = x Nf is the total number of quarks. Similarly the potential term is given as 



V 



(2vr) 



^^^(^i - k2){A+{k^))p,iA_{k,))sAdk2, 



(10) 



where V{k) is the Fourier transform of the potential V{r) given as 



V{k) = I Vif^e'^-'^dk 



(11) 



and A± are 



A^(k) = ^ (l ± -i^sin2f{k) ± a.kcos2f{k)) . (12) 



The expression for V as in Eq.( 10) can be calculated for Eq.(l) or (2a) for an effective 
potential 0. We shall however now illustrate the method with Eq.(2b) for NJL model 
corresponding to the contact potential. The total energy density then becomes 

2N r ^ ^ 

S{f} = £ = --— ^ dk\k \ cos2f - 2GN{2N + l)/^ (13) 
(2t:Y J 

with 

" sin2f{k)dk. (14) 



(27r)3 

The leading order in N here corresponds to the Hartree approximation Q]. The energy 
functional £{f} here is quadratic in f{k) and it is to be determined by minimising the 
energy density. This yields that 

tan2f{k) = ^^^^^f + ^ M/K (15) 

where, M = 2GI{2N + 1) is the dynamically generated mass. Further, substituting the 
above in Eq.( |l^) yields the self consistency relation 

2G(2A^+1) /-A M 
M = — \ — - — ^ / dk, (16) 

27r 3 J Vm^TF ^ ^ 



with A above as the ultraviolet cutoff for NJL model. Eq.(|l6]) is usually derived through 
an approximate solution to Schwinger Dyson equation 0. We followed here an alternative 
variational method with phase transition as vacuum realignment as in Eq.(5), determined 
through minimising energy density functional 0| . 

The above equation has a solution with M ^ (Goldstone phase) provided 

GA\2N + 1) > 27t\ (17) 

The energy density of | vac' > with respect to the perturbative vacuum | vac > may be 
evaluated to be 



AS = em - £ {/ = 0} _ / - y^^)* + 55j5^M^ (18) 



which is negative when a nontrivial solution to Eq.( |l^) exists or Eq.( is satisfied. The 
state with condensates | vac' > then becomes the physical vacuum. One may also calculate 
the order parameter < i/jiIj > given as 



III. GOLDSTONE THEOREM AND PION WAVE FUNCTION 

We shall now recapitulate that the present description of phase transition permits 
us to define pion also as a quark antiquark pair. From the gap equation we obtained two 
solutions for the field operators corresponding to sin2f{k) = or sin2f{k) ^ along with 
the corresponding ground state as | vac > or | vac' > respectively. 

For the case of chiral symmetry breaking, we have the gap equation 



which determines the value of the mass parameter M. Once M is determined, the function 
f{k) becomes known and hence the condensate structure of vacuum becomes known. How- 
ever, the Hamiltonian of equation (1) had chiral symmetry, which through equation (19) or 
otherwise is now seen to be broken. Hence we should have a Goldstone mode corresponding 
to a zero mass particle 0. We shall approach this theorem in a modified manner to obtain 
the wave function as a quark antiquark pair p|. When chiral symmetry remains good. 




(19) 




(20) 



Qs" I vac >= 



(21) 



where Q^"" is the chiral charge operator given as 




(22) 



For symmetry broken case however 
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Qs" I vad >^ 0. (23) 

We expect that this will describe a pion of zero total momentum. Since it will be massless it 
will also have zero energy, corresponding to the pion state. To show this we first note that 

[Q5^ ii\ = (24) 

irrespective of whether Qs" and H are written in terms of field operators corresponding to 
sin2f = or sin2f ^ since the anticommutation relation between the operators remain 
unchanged by the Bogoliubov transformation. Clearly, for the Goldstone phase, \vac' > is 
an approximate eigenstate of H with SV as the approximate eigenvalue (V being the total 
volume). With H^ff = H — SV, we then obtain from Eq.(23) that 

HeffQ^^lvac >=0 (25) 

i.e. the state Qs" | vac' > with zero momentum has also zero energy, thus corresponding to 
the massless pion. Explicitly, using Eq.(3) and Eq.(7), we then obtain with qj now as two 
component isospin doublet corresponding to (u,d) quarks above, 

I 7r"(0) >= J qi{k)\^)qi{-k)sin2f{k)dk \ vac >, (26) 

where, N^^ is a normalisation constant. The wave function for pion thus is given as propor- 
tional to u{k) = sin2f{k). The isospin and spin indices of and q for quarks have been 
suppressed. Further, with 

< tt''{6) I tt\p} >= d^'Sip), (27) 
the normalisation constant N^^ is given by 

Ar^2 ^ ^^c^ f sin^2f{k)dk = 1. (28) 



Clearly the state as in Eq.( ^6]) as the Goldstone mode will be accurate to the extent we 
determine the vacuum structure sufficiently accurately through variational or any other 
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method so that \vac' > is an eigenstate of the Hamiltonian. The above results relate pion 
wave function with vacuum structure for any example of chiral symmetry breaking. 

We may note that we could relate pion wave function to the vacuum structure since 
vacuum had an explicit structure as in equations (5). In fact, the two body condensate as 
in equations (5) for the destabilised vacuum is strictly related to the presence of the zero 
mode as seen here. 



IV. APPROXIMATE CHIRAL SYMMETRY 

While considering chiral symmetry breaking, we often use results from current algebra 
so that we may obtain numbers for approximate chiral symmetry breaking. For the sake of 
completeness, with the present mechanism, we elaborate p these results so as to use the 
same for NJL model. For this purpose we may add a small mass term to the Hamiltonian 
that breaks the chiral symmetry explicitly. Then Qs" | vac' > will not be a zero mode 
and will have finite mass. In fact the mass of the pion in the lowest order will now be m-Tr 
formally given as 

< 7r'^(0) I Hsb I vr°(0) >= m,<5(0), (29) 

where Hgi, is the symmetry breaking part of the Hamiltonian corresponding to the Hamil- 
tonian density Tigb = rmptp, m being the current quark mass. The above may be related 
to A^vr and pion decay constant as follows. Firstly we note that the identity for pion decay 
constant is 



< I I Ap) >^j^y^ . e--; (30) 
where, c-^ = 94MeV. The normalisation constant N^^ in Eq.(|26|) is then given by using 

X (5(0) = < vac' I Q%Q% \ vac' > 

= J < vac' I g% I TT^p) >dp< tt\p) I Q% I vac' >, (31) 
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where we have saturated the intermediate states with pions. The index b is summed and 
there is no summation over the index a. With Eq.(pOD and Eq.(31) we then have 



iV,-2 = l.(27r)=^ ■m^c^^ (32) 

which hnks N^^ of vacuum structure with pion mass and pion decay constant. We shall 

now substitute explicitly the pion state as in equation (26) in equation (29). We shall also 

substitute the value of normalisation constant in equation (26) by equation (32). On using 
straightforward commutation relations, we then obtain that 

< 7r'^(0) I H,, I 7r"(0) >= x -i-- < vac' \ Q\H,,Q\ \ vad > 

^ X ^Xr^ X ^ < vac I [[g%, H4 , <5%] I vac' > (33) 



m^c^2 (27r)3 2 

2 TTl — -* 

X — - < vac'lipip I vac' > x6{0). (34) 



From equation (34) and equation (29) we then obtain that 



Tfl — 

m^^ = <^i)> . (35) 

which is the familiar result for current algebra. 

V. CHARGE RADIUS OF PION 

With the wave function of the pion as above, we may next estimate the size of the 
Goldstone pion as related to the vacuum structure. The pion state with momentum p using 
translational invariance from Eq.(26) becomes 

|7r+(p) >= N^j dUi\k + ^-)\T+),f4{-k + ^-)u{k)\vad > . (36) 

In Breit frame the electric form factor is given by [|ll|] 

GE{t) = (27r)3 < 7i+i-p)\Jo\7r+{p) > (37) 

where t = —Ap^ and Jq = eip'^il). This may be evaluated directly as 



X \ui{k - p)Ui{k + p) + (A;^ - p^)u2{k - p)U2{k + p)^ . 



(38) 



In the above 



ui{k) 



\ 



{l + szn2f{\k\)) 



\k\ 



\ 



'l-mn2f{\k\)) 



(39) 



To calculate the charge radius we expand the above in powers of p and the coefficient of 
will be related to the charge radius through 

1 



GE{t) = e{l + -Rlt + ■■■). 



With GEit) as in Eq(^) we then obtain that 



- I dk 



1 2 2 



f 2 2 , 

[ k k 



where we have substituted 



(40) 



(41) 



uo{k) 



X u{k), 



(42) 



'{J sin^2f{k)dk) 
and, primes denote differentiation with respect to k. 

The above formula applies for any known vacuum realignment with condensates. Let us 
now estimate the charge radius in Nambu Jona-Lasinio model. We shall also use Eq.(34) 
for pion decay constant sothat chiral symmetry was approximately true. With Tisb = mipil), 
the extra contribution to the energy density is —m x 2NI. On extremisation the modified 
gap function is given by 

2G{2N+l)I + m M' 



tan2f{k) 



\k\ 



k 



(43) 



where M' = 2G{2N + 1)/ + m satisfies the equation parallel to Eq.(|T^) given as 

dk 



2G{2N+1) M 
(2^0^ ^ 



10 



+ m/M'. 



(44) 



Thus here we may also have a vacuum reahgnment. 

We shall now choose an optimal set of parameters for GA^, A and m. Then, for example, 
with A = 420 MeV, GA2=2.24 and m=16 MeV, we get M=305 MeV, < -tptp >^/^= 220 
MeV, i?^=0.25 fermi^ and 0,^= 94 MeV. Here we have taken = 138 MeV. As a further 
illustration to see how the the corresponding quantities change with parameters, for A=500 
MeV, GA2=2.15, m=10 MeV, we have, M=320 MeV, < -^^^ >i/3=255 MeV, Rl=0.20 
fermi^ and Ct^=93 MeV. We note that the pion structure as arising from vacuum realignment 
appears to give a smaller value of charge radius than would be expected. In fact, with GA^ 
= 2.0, A = 700 MeV and m= 5 MeV g, we obtain that M= 360 MeV, < -tjj^ >y^= 342 
MeV, R^= 0.13 fermi^ and C7r=103 MeV. Thus the above set of parameters do not appear 
to be acceptable [|16|. We may also note that the four component Dirac field operators 



for the quarks will change the above numbers as examined elsewhere |]T2[, which however 
does not change the above remarks. The above illustrates the nature of constraints derived 
for symmetry breaking through determination of pion wave function. A parallel approach 
with Bogoliubov transformations and Schwinger Dyson equation has been used to obtain 
Salpeter wave function for the pion, which, however does not permit the definition of pion 
as a state since the wave function is not normalisable and therefore did not give rise to such 
constraints. 



VI. NEW MODES IN VACUUM 

When vacuum has a structure, there can be excitations present due to such a structure. 
For chiral symmetry breaking, let us substitute 

^^' + Mla{x), (45) 

where Msc is a mass parameter and (y{x) may represent the scalar field of vacuum fluctu- 
ations. Then a{x) can represent quantum fluctuations of the condensate. In fact, we may 
evaluate 
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< vac'\{tp{x)'ilj{x) - ^j!^){ip{y)ij{y) - /i^)|mc' > 
Ml < vad\a{x)a(v]\vad > 



Mt 



Jk.(x-y) 



-dk. 



{271 f J 2(A;2 + m2)i/2 
where we approximate <y{x) by a free field of mass mg-. Let us define 

I{k) = / dxexp{-ik.x) < vac'\{ijj{x)'ip{x) - /i^)(-?/'(O)'0(O) - fj,^)\vac' > 



In that case, clearly free field approximation for a{x) corresponds to 

I{k) = , . 
2\lml + P 

Explicit evaluation of the left hand side of equation 



(46) 



(47) 



(48) 



in the limit of small \k\ yields 



I{k) ^ - 



3 ^M' 



-Pi 



A' 



+ 



+ 



AM^ 



_ 3 -1. A , 

l6(A2 + M2)2 ' 8'(A2 + M2)2 ' 16'(A2 + M2) ~ 16' ^^"^ ^M' 



1 



in-" 



Mi Mi 
- k.- 



(49) 



.2m^ 4m^. 

where, in equation (^) we have kept terms upto k"^. Equating equal powers of k in equation 
(^) and eliminating Mgc in favor of rricr yields, with x = M/A, 
o A' 



% — + xHan H-) 

3 ^ x' 



1 1 



lb 



I 3 I l_xP__ 

x' ^ l& (l+x2) 8 (l+x2)2 



(50) 

. (l+a;2)2 

We may next estimate the mass of such a mode for different values of A and M obtained 
in the previous section. For example, for A = 420 MeV and M = 305 MeV, vrL„ ^ 2.07M; 
for A = 500 MeV and M = 320 MeV, ~ 2.27M and for A = 700 MeV and M = 360 
MeV, nif, ~ 2.65M. These may be compared with the mass of a field obtained through an 
approximate determination of the pole of the propagator with polarisation insertion, which 
is given as = (4M^ + m^)^/^ M. 



VII. DISCUSSIONS 



We thus consider here chiral symmetry breaking as a vacuum realignment with an explicit 
construct for destabilised vacuum. The new feature of this approach ^ is that it enables 
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us to relate the function that describes the vacuum structure determined variationally to 
the wave function of the pion as the locahsed Goldstone mode in a straightforward manner. 
This language is not only physically appealing reproducing the conventional results but also 
puts severe constraints on the parameters for symmetry breaking as illustrated here for NJL 
model. Some other aspects of low energy hadronic properties as related to the vacuum 
structure for chiral symmetry breaking have been discussed elsewhere [O . 



ACKNOWLEDGMENTS 

The authors are thankful to Snigdha Mishra, Amruta Mishra, S. N. Nayak and P. K. 
Panda for many useful discussions. SPM acknowledges to Department of Science and Tech- 
nology, Government of India for the research grant SP/S2/K-45/89 for financial assistance. 



13 



REFERENCES 



[1] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122, 345 (1961); 124, 246 (1961). 

[2] J.R. Finger and J.E. Mandula, Nucl.Phys.B199, 168 (1982); S.L. Adler and A.C. Davis, 
Nucl. Phys. B244, 469 (1984); R. Alkofer and P. A. Amundsen, Nucl. Phys.B306, 305 
(1988). 

[3] S. Li, R. S. Bhalerao and R. K. Bhaduri, Int. J. Mod. Phys. A6, 501 (1991), V. Bernard, 
Phys. Rev. D34, 1601 (1986). 

[4] J. da Providencia, M. C. Ruivo and C. A. de Sousa, Phys. Rev. D36, 1882 (1987); B. 
Haeri and M. B. Haeri, Phys. Rev. D43, 3732 (1991); B. Rosenstein, B.J. Warr, and 
S.H. Park, Phys. Rep. 205, 59 (1991); U. Vogl and W. Weise Prog. Part. Nucl. Phys. 
27, 195 (1991); L. N. Chang, N. P. Chang and K. C. Chou, Phys Rev D43, 596 (1991). 
Lay Nam Chang and Ngee Pong Chang, Phys. Rev. D45, 2988 (1992); S.P. Klevensky, 
Rev. Mod. Phys. 64, 649 (1992). 

[5] H. Mishra, S.P. Misra and A. Mishra, Int. J. Mod. Phys. A3, 2331 (1988); A. Mishra, 
H. Mishra, S. P. Misra and S. N. Nayak, Phys. Lett.251B, 541 (1990); ibid, Phys. Rev. 
D44, 110 (1991); ibid, Pramana(J. Phys.)37, 59 (1991) ibid, Z. Phys. C57, 233 (1993). 

[6] A. Mishra, H. Mishra and S. P. Misra, Z. Phys. C57, 241 (1993). 

[7] J. Goldstone, Nuov. Cim. 19, 154 (1961); J. Goldstone, A. Salam and S. Weinberg, 
Phys. Rev. 127, 965 (1962). 

[8] S. P. Misra, Talk given at Symposium on Quantum Field Theory and Statistical Me- 
chanics, Calcutta, 1992, |hep-ph/9212293| . 



[9] Y. Nambu in New Theories in Physics, Proceedings of the Eleventh International Sym- 
posium on Elementary Particle Physics, Kazimierz, Poland, edited by Z. Ajduk, S. 
Pokoroski and A. Trautman (World Scientific, Singapore, 1989), W. J. Marciano, Phys. 



14 



Rev. Lett. 62, 2793 (1989); Phys. Rev. D41, 219 (1990), W. A. Bardeen, C. T. Hill and 
M. Linder, Phys. Rev. D41, 1647 (1990). 

[10] K.S. Babu and Rabindra N. Mohapatra, Phys. Rev. Lett. 66, 556 (1991) 

[11] S. P. Misra, Phys. Rev. D18, 1661 (1978); ibid, D18, 1673 (1978). 

[12] A. Mishra and S. P. Misra, Z. Phys. C 58, 325 (1993). 

[13] A. Le Yaouanc, L. Ohver, S. Ono, O. Pene and J.C. Raynal, Phys. Rev. Lett. 54, 506 
(1985). 

[14] H. Mishra, S. P. Misra, P. K. Panda and B. K. Parida, Int. J. Mod. Phys. El, 405 
(1992). 

[15] J.J. Sakurai in "Currents and Mesons", The university of Chicago press 1969 p. 83 

[16] There is a disagreement over the value of A. Further, Eq.(17) and the corresponding 
equation by Bernard in Ref.3 have a difference of a factor 2. We presume that this is 
a printing error since GA^ they have taken does not satisfy their bound where as it is 
consistent with Eq.(17). 

[17] A. Amer, A. Le Yaouanc, L. Oliver, O. Pene and J.C. Raynal, Phys. Rev. Lett. 50, 87 
(1983); ibid, Phys. Rev. D28, 1530 (1983). 



15 



